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Abstract 

We review on the viability of _F(_R)-gravity. We show that recent cosmic acceleration, 
radiation/matter-dominated epoch and inflation could be realized in the framework 
of _F(i?)-gravity in the unified way. For some classes of -F(_R)-gravity, the correction 
to the Newton law is extremely small and there is no so-called matter instability (the 
very heavy positive mass for additional scalar degree of freedom is generated). The 
reconstruction program in modified gravity is also reviewed and it is demonstrated 
that any time-evolution of the universe expansion could be realized in F(i?)-gravity. 
Special attention is paid to modified gravity which unifies inflation with cosmic ac- 
celeration and passes local tests. It turns out that such a theory may describe also 
dark matter. 



1 Introduction 

Recent astrophysical observations indicate that the accelerating expansion of the universe has started 
about five billion years ago and the present universe is flat. This implies the existence of dark energy, 
that is, unknown component in the universe. 

Usually the evolution of the universe can be described by the FRW equation: 
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H 2 = P . (1) 



Here the spatial part of the universe is assumed to be flat. We denote the Hubble rate by H, which is 
defined in terms of the scale factor a by 

In (JTJ) , p expresses the energy density of the usual matter, dark matter, and dark energy. The dark energy 
could be cosmological constant and/or a matter with 'equation of state (EoS)' parameter w, which is less 
than —1/3 and is defined by 

u; = £<-l/3. (3) 
P 

Instead of including unknown exotic matter or energy, one may consider the modification of gravity, 
which corresponds to the change of the l.h.s. in (JTJ) . 

An example of such modified gravity pretending to describe dark energy could be the scalar-Einstein- 
Gauss-Bonnet gravity [T], whose action is given by 

1 „ 1 



Here Q is Gauss-Bonnet invariant: 

Q = R 2 - AR^R^ + R^pvR^P* . (5) 



1 E-mail:nojiri@phys. nagoya-u.ac.jp 

2 E-mail:odintsov@ieec. uab.es, also at Lab. Fundam. Study, Tomsk State Pedagogical University, Tomsk 



Another example is so-called F(i?)-gravity (for a review, see 0). In F(i?)-gravity models, the scalar 
curvature R in the Einstein-Hilbert action 



Seh = / d xy^gR , (6) 

is replaced by a proper function of the scalar curvature: 

S F (R) = J d 4 x^F(R) . (7) 

Recently, an interesting realistic theory has been proposed in [3J, where F(R) is given by 

1 m 2 c (R/m 2 ) n 

F(R) — 7T^ {R + fHs(R)) , f HS (R) = 7^ ■ (8) 

In this model, R is large even in the present universe, and fns(R) could be expanded by the inverse 
power series of R: 

m 2 Cl m 2 Cl / R 



MR)~-^ + ^{- 2 ) + (9) 

Then there appears an effective cosmological constant A c ff as A c fr = m 2 ci/c 2 , which generates the 
accelerating expansion in the present universe 

In the HS-model ([8]), there occurs a flat spacetime solution, where R = 0, since the following condition 
is satisfied: 

Bmf H s(R)=0. (10) 
a — >u 

An interesting point in the HS model is that several cosmological conditions could be satisfied. 

In the next section, we review on the general properties of i 7 '(i?)-gravity. After some versions of 
-F(i?)-gravity were proposed as a model of the dark energy, there appeared several criticisms/viability 
criteria, which we review in Section [3J It is shown how the critique of modified gravity may be removed 
for realistic models. In Section [51 we propose models [3] and [5], which unify the early-time inflation 
and the recent cosmic acceleration and pass several cosmological constraints. Reconstruction program 
for i 7 '(i?)-gravity is reviewed in Section [5j The partial reconstruction scenario is proposed. Section six 
is devoted to the description of dark matter in terms of viable modified gravity where composite scalar 
particle from F(R) gravity plays the role of dark particle. Some summary and outlook is given in the 
last section. 



2 General properties of -gravity 



In this section, the general properties of the F(i?)-gravity are reviewed. For general F(i?)-gravity, one can 
define an effective equation of state (EoS) parameter. The FRW equations in Einstein gravity coupled 
with perfect fluid are: 

p=^H 2 , p = ~ (ZH 2 + 2H) . (II) 
For modified gravities, one may define an effective EoS parameter as follows: 

WcS = 1 ff^ • (12) 

The equation of motion for modified gravity is given by 

\g^F(R) - R pu F'(R) - g^UF'{R) + V^V U F'(R) = -yT (m))1 „ . (13) 
By assuming spatially flat FRW universe, 

ds 2 = -dt 2 + a{tf ( dxl ) 2 > ( 14 ) 

i=l,2,3 
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the FRW-like equation follows: 



= + 3 (h 2 + H) F'(R) - 18 Uh 2 H + HH) F"{R) + K 2 p (m) (15) 



There may be several (often exact) solutions of (fT5|) . Without any matter, assuming that the Ricci 
tensor could be covariantly constant, that is, i? M „ oc g^ v , Eq. (fT3|) reduces to the algebraic equation: 

= F{R) - 2RF{R) . (16) 

If Eq. (|16[) has a solution, the Schwarzschild (or Kerr) - (anti-)de Sitter space is an exact vacuum solution 
(see[6] and refs. therein). 

When F(R) behaves as F(R) oc R m and there is no matter, there appears the following solution: 

(m-l)(2m-l) 

H 2p2 , (17) 

which gives the following effective EoS parameter: 

6m 2 — 7m — 1 
3(m — l)(2m — 1) 

When F(R) cx R m again but if the matter with a constant EoS parameter w is included, one may get 
the following solution: 

2m 



H~*&±V, (19) 



and the effective EoS parameter is given by 



S 



w + 1 , s 

w eS = -l + . (20) 

m 

This shows that modified gravity may describe early/late-time universe acceleration. 

3 Criticism of -gravity 

Just after the -F(i?)-models were proposed as models of the dark energy, there appeared several works 
[3 [8] (and more recently in [21 [10]) criticizing such theories. 

First of all, we comment on the claim in [7]. Note that one can rewrite F(i?)-gravity in the scalar- 
tensor form. By introducing the auxiliary field A, we rewrite the action ([7} of the i ? (i?)-gravity in the 
following form: 

1 J d 4 xV=d {F'(A) {R-A) + F(A)} . (21) 

By the variation over A, one obtains A — R. Substituting A = R into the action (|21|). one can reproduce 
the action in |7|). Furthermore, we rescale the metric in the following way (conformal transformation): 

g^^e^g^, a = -]nF'(A), (22) 

Hence, the Einstein frame action is obtained: 

Se = ^ J d 4 x^^R~^d p ad a a-V(a^j , 

V{a) = eM^)-e-/(H^)H^-^ (23) 

Here g (e _CT ) is given by solving the equation a = — In (1 + f'{A)) = lnF'(^4) as A = g (e _cr ). Due to the 
scale transformation ([22]) . there appears a coupling of the scalar field a with usual matter. The mass of 
a is given by 

a _ 1 d»yfr) _ 1 f A 4F(A) 1 ] 

a ~2 da 2 2\F'(A) (F'(A)) 2 F"{A) J ' 1 ' 
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Unless m a is very large, there appears the large correction to the Newton law. Naively, one expects the 
order of the mass m a could be that of the Hubble rate, that is, m a ~ H ~ 10 -33 eV, which is very light 
and the correction could be very large, which is the claim in [7J. 

We should note, however, that the mass m a depends on the detailed form of F(R) in general [IT] . 
Moreover, the mass m a depends on the curvature. The curvature on the earth -R ear th is much larger than 
the average curvature i? so iar in the solar system and i£ S olar is also much larger than the average curvature 
in the unverse, whose order is given by the square of the Hubble rate H 2 , that is, i? e arth ^ -Rsoiar ^ H 2 . 
Then if the mass becomes large when the curvature is large, the correction to the Newton law could be 
small. Such a mechanism is called the Chameleon mechanism and proposed for the scalar-tensor theory 
in [12]. In fact, the HS model [3] has this property and the correction to the Newton law can be very 
small on the earth or in the solar system. In the HS model, the mass m a is given by (see also |13j ) 

TO 2 C 2 / R \n+2 

™* ~ 2n(n+\) Cl [m 2 ) ' (25) 

Here the order of the mass-dimensional parameter m 2 could be m 2 ~ 10 _64 eV 2 . Then in solar system, 
where R ~ l(T 61 eV 2 , the mass is given by ml ~ 10~ 58+3 ™eV 2 and in the air on the earth, where 
R ~ l(T 50 eV 2 , m 2 a ~ icr 36 + 14 " eV 2 . The order of the radius of the earth is 10 7 m - (lCT^cV) -1 . 
Therefore the scalar field a could be heavy enough if n 3> 1 and the correction to the Newton law is not 
observed being extremely small. On the other hand, in the air on the earth, if we choose n = 10, for 
example, one gets the mass is extremely large: 

m a ~ 10 43 GeV ~ 10 29 x M P i anck . (26) 

Here Mpianck is the Planck mass. Hence, the Newton law correction should be extremely small. 

Let us discuss the matter instability proposed in [5], which may appear when the energy density or 
the curvature is large compared with the average one in the universe, as is the case inside of the planet. 
Multiplying g^ v with Eq. (fT3"]) , one obtains 



F^{R) F'(R)R 2F(R) n. 

M+ fP)(fl) V ^ V R+ 3F^)'W^)-6fTW) T - (27) 



Here T is the trace of the matter energy-momentum tensor: T = T,j\ . We also denote d n F(R)/dR n 

by FW(R), Let us now consider the perturbation from the solution of the Einstein gravity. We denote 
the scalar curvature solution given by the matter density in the Einstein gravity by ~ (k 2 /2)p > 
and separate the scalar curvature R into the sum of i?b and the perturbed part R p as R = R(, + R p 
(\R P \ <§; \Rb\)- Then Eq. (|2"T]) leads to the perturbed equation: 

° DRb+ FT^r) VpRbV Rb+ 3F( 2 )(R b ) 

2F[Rb) ' h 1 DR P + 2^^-W p R b W"R p + U(R b )R p . (28) 



3F( 2 )(i? 6 ) 3FW(Rb) p FW(R, 
Here U(R b ) is given by 



_ ( FWjRt) F^(R b ) 2 \ R b 
= {F^iRV) ~ F<W ) P " + "3" 

F^(R b )F^(R b )R b F^(R b ) , 2F(R b )F^(R b ) F&(R b )R b 



3F( 2 )(R b ) 2 3F( 2 ){R b ) 3F( 2 ){R b ) 2 3F( 2 )(R b ] 
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(29) 



It is convinient to consider the case that R b and R p are uniform, that is, they do not depend on the 
spatial coordinate. Hence, the d'Alembertian can be replaced with the second derivative with respect to 
the time coordinate: OR p — > —d 2 R p and Eq. (l29]) has the following structure: 

= -d 2 R p + U(R b )R p + const. . (30) 
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Then if U (Rb) > 0, R p becomes exponentially large with time t: R p ~ eV^W' anc [ the system becomes 
unstable. In the l/i?-model |14j , since the order of mass parameter m M is 

~ 10 18 sec ~ (lO- 33 eV) _1 , (31) 

one finds 

[/( J R h )^- J R h + #|^^^(lO- 26 sec) 



6/i 4 fj, 4 \gcm~ 3 



i? b ^(l0 3 scc)- 2 f^4) (821 



Hence, the model is unstable and it would decay in 10 26 sec (for planet size). On the other hand, in 
1/R+ i? 2 -model [TQ, we find 

U(Ro) ~ ^- > . (33) 

Then the system could be unstable again but the decay time is ~ 1, 000 sec, that is, macroscopic. In HS 
model |3i . U{Rb) is negative [15]: 

U(R Q )^-^±^f<0. (34) 
c\n\n + 1) 

Therefore, there is no matter instability [13] . 

Let us discuss the critical claim against modified gravity in [9] [TP] . As shown in ([T6]) . as an exact 
solution, there appears de Sitter-Schwarzschild spacetime in F(i?)-gravity. The claim in [9] [10] is that the 
solution does not match onto the stellar interior solution. Since it is difficult to construct explicit solution 
describing the stellar configuration even in the Einstein gravity, we now proceed in the following way: 
First, we separate F(R) into the sum of the Einstein-Hilbert part and other part as F(R) = R + f(R). 
Then Eq. ([T3]) has the following form: 

= -\g»u (f(R) + A) + R^f'(R) + g^af(R) - V^V„/'(iZ) . (35) 

Here —A is the value of f(R) in the present universe, that is, A is the effective cosmological constant: 
A = —f(Ro). We now treat the r.h.s. in ([55]) as a perturbation. Then the last two derivative terms in 
(l35]l could be dangerous since there could be jump in the value of the scalar curvature R on the surface of 
stellar configuration. Of course, the density on the surface could change in a finite width A as in Figure 
Q] and the derivatives should be finite and the magnitude could be given by 

d„ ~ i ■ (36) 

One now assumes the order of the derivative could be the order of the Compton length of proton: 

d p ~ m p ~ 1 GeV - 10 9 eV (37) 

Here m p is the mass of proton. It is also assumed 

R e - 10~ 47 eV 2 , (38) 

that is, the order of the scalar curvature R is given by the order of it inside the earth. 
In case of the 1 /R model [2] , one gets 

2 4 

in Uj 

Uf'{R) ~ \7^ v f'(R) ~ -g- ~ 10- 20 eV 2 » R e . (39) 

Then the perturbative part could be much larger than unperturbative part in (|35p . say, R ~ R e ~ 
10 _47 eV 2 . Therefore, the perturbative expansion could be inconsistent. 
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Figure 1: Typical behavior of R and p near the surface of the stellar configuration. 
In case of the model [3], however, we find 

Of'(R) ~ VpVvf'iR) ~ ^ ( 4 ^ 10 " 3 " 17 " eV ' ■ (4°) 



Then if n > 2, we find □/'(i?), V v f (R) <C R e and therefore the perturbative expansion could 
be consistent. This indicates that such modified gravity model may pass the above test. Thus, it is 
demonstrated that some versions of modified gravity may easily pass above tests. 

4 Unifying inflation and late-time acceleration 

In this section, we consider an extension of the HS model [5] to unify the early-time inflation and late-time 
acceleration, following proposals [H [5] . 

In order to construct such models, we impose the following conditions: 



Condition that inflation occurs: 



lim f(R) = -A, . (41) 

H — >oo 



Here Aj is an effective early-time cosmological constant. 
Instead of (|4ip one may impose the following condition 

lim f(R) = aRT . (42) 

R — >OG 

Here m and a are positive constants. Then as shown in (fi"9)) . the scale factor a(t) evolves as 

(43) 

and w c ff = — 1 + 2/3/iq. Here w is the matter EoS parameter, which could correspond to dust or 
radiation. We assume m ^S> 1 so that H/H 2 ^> 1. 

• The condition that there is flat spacetime solution is given as 

/(0) = (44) 

• The condition that late-time acceleration occurs should be 

f(R ) = -2Ro , f'(Ro) ~ . (45) 

Here i?o is the current curvature of the universe and we assume Rq > Ro- Due to the condition 
(|45)l . f(R) becomes almost constant in the present universe and plays the role of the effective small 
cosmological constant: A; ~ —f(Ro) = 2Rq. 



(» 



f(R) 




Figure 2: The typical behavior of f(R) which satisfies the conditions (l4Tj) . (|44|). and (|45ll . 
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Figure 3: The typical behavior of /(-R) which satisfies the conditions l)42p. (|44ll . and (j45 




Figure 4: The potential in the inflationary epoch. 



The typical behavior of f(R) which satisfies the conditions (|4"Tj) . (T4"4"|) , and (|4"5| is given in Figure [5] 
and the behavior of f(R) satisfying (|4Tj) . (|42|) . and (|45|) is given in Figure [TJ 

Some examples may be of interest. An example which satisfies the conditions (j4Tj) . (|44|) . and (|45|) is 
given by the following action [4]: 

/ p R x2ri+l , R 2n+1 

/(i?) = ^=-^ ±^ T • (46) 

h + h{(R-Ro) 2n+1 + R 2 o n+1 } 

Here n is a positive integer. The conditions (|4^]) and ([43]) require 

n2n+l i 

-2^o, - r =A i . (47) 



One can now investigate how the exit from the inflation could be realized in the model (|46p . It is easier 
to consider this problem in the scalar-tensor form (Einstein frame) in (|23p . In the inflationary epoch, 
when the curvature R — A is large, f(R) has the following form: 



Hence, one gets 



and 



( 2 " +1 )/° (49 ) 



2n + l 



1_ _ 2(n+l)/ (_fi°_\ 2 " +2 
A /i 2 U2n + l)/ , 



Note that the scalar field a is dimensionless now. Let us check the condition for the slow roll, \V /V\ <C 1. 
Since 

V(a) \(2n+l)f J 

if we start with a ~ 1, one finds 

which is very small and the slow roll condition is satisfied. 

Thus, the value of the scalar field a increases very slowly as in Figure @] and the scalar curvature R 
becomes smaller. When er becomes large enough and R becomes small enough, the inflation could stop. 



Another possibility to achieve the exit from the inflation is to add small non-local term to gravitational 
action. 

We now consider another example, where f(R) satisfies the conditions (|4"2"j) . l|34")). and (|4"S"jl [5]: 

Here a, /?, and 7 are positive constants and n is a positive integer. When the curvature is large (R — ► 00), 
/(i?) behaves as 

f{R) - -i?" ■ (54) 
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To achieve the exit from the inflation, more terms could be added in the action. Since the derivative of 
f(R) is given by 

(/ , ni?™" 1 (ajR 2n - 2aR n - (3) , 

f (R) = — 5 — , 55) 

(l + 7# n ) 

we find the curvature Rq in the present universe, which satisfies the condition f'(Ro) — 0, is given by 

^{Ul + J^S)]^ , (56) 



7 V V a 



and 



7 2 V 2 + + /? 7 /a / 

Let us check if we can choose parameters to reproduce realistic cosmological evolution. As a working 
hypothesis, we assume (3"f/a 3> 1, then 

Ro ~ (/?/« 7 ) 1/2n , /(i? ) = -2i?o ~ -y9/T (58) 

We also assume f(Ri) ~ (a/j)R™ ~ i?/. Here i?j is the curvature in the inflationary epoch. As a result, 
one obtains 

a ~ 2i? i?Q 2 ", /3 - ARlR^ 2n R^-\ 7 - 2R R^ 2n R^ 1 . (59) 
Hence, we can confirm the assumption (3j/a 1 if n > f as 

^ ~ 4i?gi?Q 2 "i?f- 2 ~ lO 228 *"- 1 ) » f . (60) 

Thus, we presented modified gravity models which unify early-time inflation and late-time acceleration. 
One should stress that the above models ([46]) and (|53[) satisfy the cosmological constraints/local tests in 
the same way as in the HS model [3]. 

5 Reconstruction of gravity 

In this section, it is shown how we can construct F(R) model realizing any given cosmology (including 
inflation, matter-dominated epoch, etc) using technique of ref . [15] . The general i ? (i?)-gravity action with 
general matter is given as: 

S = J d A xV=jj {F(R) + £ matter } . (61) 
The action (|61|) can be rewritten by using proper functions P(4>) and Q{4>) of a scalar field 4>: 

S = J d^X^—g {P(4>)R + Q{4>) + ^matter} . (62) 

Since the scalar field <f> has no kinetic term, one may regard <fi as an auxiliary scalar field. By the variation 
over </>, we obtain 

Q = P'^)R + Q'{4>) , (63) 
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which could be solved with respect to <f> as 4> — 4>{R)- By substituting cj> — 4>(R) into the action (JHU), we 
obtain the action of -F(i?)-gravity where 

F(R)=P(<p(R))R + Q(<f>(R)) . (64) 

By the variation of the action (|62p with respect to g^ v , the equation of motion follows: 

o = -\g»» {P{<t>)R + QW} - R„ V P{4>) + v M v„p(0) - 5 ^v 2 p(0) + ir^ (65) 

In FRW universe ((HJ), Eq. ([63)l has the following form: 

= -6ff 2 P(0) - Q{4>) - 6H dP{ f} t)} + p 

at 

= U H + 6HAp W + Q^ + 2 d ^^ + 4H d -^l +P (66) 
V / at at 

By combining the two equations in (|66[) and deleting Q(<j>), we obtain 

d 2 PU(t)) dP(6(t)) ■ , , s 

= 2 — tpt 11 - 211 — T^+^+P + r 67 

ar^ at 

Since one can redefine properly as </> = (j>(<p), we may choose </> to be a time coordinate: cj> — t. Then 
assuming p, p could be given by the corresponding sum of matter with a constant EoS parameters w, 
and writing the scale factor a(t) as a = aoe 9 '*' (a : constant), we obtain the second rank differential 
equation: 

= " 2 ^)^T + 4 9"WPW + £ C 1 + "0 ftoao 3(1+t " <) e- 3(1+tUi)s(<A) • (68) 

If one can solve Eq. (|6"8)) . with respect to P(<f>), one can also find the form of Q(4>) by using (f6"6")) as 

Q(<P) = -6 (, 9 '(0)) 2 P(0) - 6 5 '(0)^ + £ P,o«o 3(1+tOl) e- 3(1+tUl)s(0) • (69) 

Thus, it follows that any given cosmology can be realized by some specific P(i?)-gravity. 
We now consider the cases that (|68[) can be solved exactly. A first example is given by 

g'{4>) = go + § . (70) 

For simplicity, we neglect the contribution from matter. Then Eq. (|68|) gives 

d 2 P {_ gA dP 2gt 

~§)~d4~~$ 



The solution of ([711) is given in terms of the Kummer functions or confluent hypergeometric functions: 

P = z a F K (a, r , z) , z x -^F K {u - 7 + 1, 2 - 75 z) (72) 

Here 



1 + 5i ± v/.g 2 + 10.9i + 1 
z = g (j> , a= 



1± yM±MH t F K (a, r> z) = t a{ r + 1 l'"l a + n ~ 1) A 



(73) 



Eq. ([70| gives the following Hubble rate: 



H = .90 + f ■ (74) 
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Then when t is small, H behaves as 

(75) 

which corresponds to the universe with matter whose EoS parameter is given by 

u, = -l+A (76) 

On the other hand, when t is large, we find 

H -> go , (77) 

that is, the universe is asymptotically deSitter space. 

We now show how we could reconstruct a model unifying the early-time inflation with late-time 
acceleration. In principle, one may consider g(4>) satisfying the following conditions: 

• The condition for the inflation (t = <f> — ► 0): 

g"(0) = , (78) 

which shows that H(0) = g'(0) is almost constant, which corresponds to the asymptotically deSitter 
space. 

• The condition fot the late-time acceleration (at t = <j> ~ to): 

g"(t ) = o, (79) 

which corresponds to the asymptotically deSitter space again. 
An example could be 

g / W=go + gi \" ^ " ■ (80) 

Here g^ 1 g%, and c are positive constants and n is positive integer greater than 1. Note that g (<f>) is a 
monotonically decreasing function of <fr if < < t We also assume 

< g - « ,9o ■ (81) 

c 

One should note that g'(0) — go corresponds to the large Hubble rate in the inflationary epoch and 

g'ito) = go — 91 Q to the small Hubble rate in the present universe. It is very difficult to solve (j68| with 
(|80|) . so we expand g'{4>) for small <j>. For simplicity, we consider the case that n = 2 and no matter 
presents. Then 

g(cb)= 9 0~^<i> 2 + O(^orgl) . (82) 



4 + c- 



Hence, one gets 



P{4>) = Po + Pie 90 ^- 

l e ffo<t> _ p 3 

9o 



2 9l t 



3 50 go So J I go go 



4 + c 

+ 0(g() . (n:-!) 



Using boundary conditions we can specify different modified gravities which unify the early-time inflation 
with late-time acceleration. The important element of above reconstruction scheme is that it may be 
applied partially. For instance, one can start from the known model which passes local tests and describes 
the late-time acceleration. After that, the reconstruction method may be applied only at very small times 
(inflationary universe) to modify such a theory partially. As a result, we get the modified gravity with 
necessary early-time behavior and (or) vice-versa. 
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6 Dark Matter from F(i?)-gravity 



It is extremely interesting that dark matter could be explained in the framework of viable .F(i?)-gravity 
which was discussed in previous sections. 

The previous considerations of F(i?)-gravity suggest that it may play the role of gravitational alter- 
native for dark energy. However, one can study -F(i?)-gravity as a model for dark matter. There have 
been proposed several scenarios to explain dark matter in the framework of F(i?)-gravity. In most of such 
approaches [IB], the MOND-like scenario or power-law gravity have been considered. In such scenarios, 
the field equations have been solved and the large-scale correction to the Newton law has been found and 
used as a source of dark matter. 

There was, however, an observation [17] that the distribution of the matter is different from that of 
dark matter in a galaxy cluster. From this it has been believed that the dark matter can not be explained 
by the modification of the Newton law but dark matter should represent some (particles) matter. 

It is known that -F(i?)-gravity contains a particle mode called 'scalaron', which explicitly appears when 
one rewrites -F(i?)-gravity in the the scalar-tensor form (|23j> . In the Einstein gravity, when we quantize 
the fluctuations over the background metric, we obtain graviton, which is massless tensor particle. In 
case of F(.R)-gravity, when one quantizes the fluctuations of the scalar field in the background metric, one 
gets the massive scalar particles in the addition to the graviton. Since the scalar particles in i^(i?)-gravity 
are massive, the pressure could be negligible and the strength of the interaction between such the scalar 
particles and usual matter should be that of the gravitational interaction order and therefore very small. 
Hence, such scalar particle could be a natural candidate for dark matter. 

In the model [3] or our models (j46|) and ([53]) . the mass of the effective scalar field depends on the 
curvature or energy density, in accord with so-called Chameleon mechanism. As our models (I46[) and 
(|53|) describe the early-time inflation as well as late-time acceleration, the 'scalaron' particles, that is, 
the scalar particles in i 7 '(i?)-gravity, could be generated during the inflationary era. An interesting point 
is that the mass could change after the inflation due to Chameleon mechanism. Especially in the model 
(|46|) . the mass decreases when the scalar curvature increases as shown in (|49|) . Hence, in the inflationary 
era, when the curvature is large, one may consider the model where m a is large. After the inflationary 
epoch, the scalar particles, generated by the inflation, could lose their mass. Since the mass corresponds 
to the energy, the difference between the mass in the inflationary epoch and that after the inflation could 
be radiated as energy and could be converted into the matter and the radiation. This indicates that the 
reheating could be naturally realized in such model. Let the mass of a in the inflationary epoch be mi 
and that after inflation be m A - Then for N particles, the radiated energy En may be estimated as 

E=(m I -m A )N, (84) 

which could be converted into radiation, baryons and anti-baryons (and leptons). It is believed that the 
number of early-time baryons and anti-baryons is 10 10 times of the number of baryons in the present 
universe. Since the density of the dark matter is almost five times of the density of the baryonic matter, 
we find 

mi > l0 10 m A . (85) 

In the solar system, one gets A = R ~ 1CT 61 eV 2 . Then if n > 10 ~ 12 and A., ~ io 20 ~ 38 , the order 
of the mass m a is given by 

m 2 _ 1() 239~295-10n eV 2 ^ (86) 

which is large enough so that a could be Cold (non-relativistic) Dark Matter. On the other hand, in 
l/i?-model, the corresponding mass is given by 

m^-^-lO-^eV 2 . (87) 

Here fi is the parameter with dimension of mass and /i ~ 10 _33 eV. The mass m\/n is very small 
and cannot be a Cold Dark Matter. The corresponding composite particles can be a Hot (relativistic) 
Dark Matter but Hot Dark Matter has been excluded due to difficulty to generate the universe structure 
formation. 
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In the inflationary era, the spacetime is approximated by the de Sitter space: 

£=1,2,3 

Then the scalar particle a could be Fourier-transformed as 



ds 2 = -dt 2 + e 2H ° t J2 (^T ■ ( 88 ) 



a = J d 3 ka(k,t)e- lk x . (89) 
Hence, the number of the particles with k created during inflation is proportional to e v ~* . Here 




(90) 



Then if 

ml 9 



H 2 > 4 ' < 91 > 

sufficient number of the particles could be created. 

In the original /(i?)-frame (J7]), the scalar field a appears as composite state. The equation of motion 
in /(i?)-gravity contains fourth derivatives, which means the existence of the extra particle mode or 
composite state. In fact, the trace part of the equation of motion (|13|) has the following Klein-Gordon 
equation-like form: 

3V 2 f'(R) = R + 2f(R)-Rf'{R)-K 2 T. (92) 

The above trace equation can be interpreted as an equation of motion for the non trivial 'scalaron' f'(R). 
This means that the curvature itself propagates. In fact the scalar field a in the scalar-tensor form of the 
theory can be given by 'scalaron', which is the combination of the scalar curvature in the original frame: 

a = -ln(l + f'(R)) . (93) 

Note that the 'scalaron' is different mode from graviton, which is massless and tensor. 

Eq. (l4T))) shows that the mass, which depends on the value of the scalar field a, is given by 



m 2 ~ k ( 2 J^±1) ( ft \ ^ a -m (94) 
f?\2n + 2j {(2n + l)f > 



If the curvature becomes small, a becomes large and m 2 a decreases. Then the scalar particles lose their 
masses after the inflation. The difference of the mass in the inflationary epoch and that after the inflation 
could be radiated as energy and can be converted into the matter and the radiation. 
By substituting the expression of <r |49|) into ([94]) . one obtains 



ml ~ 2(2n+l)(n + l)/ ' (95) 

Note that A corresponds to the scalar curvature. Let denote the value of A in the inflationary epoch by 
A j and that after the inflation by A a- Then the condition ([55)) shows 



mj ( Ai 



m A \A A 



n+3/2 



> 10 iU . (96) 



For the model with n = 2, the condition ([85| or (|96|) could be satisfied if Aj/Aa > 10 3 , which seems to 
indicate that the reheating could be easily realized in such a model. 

Now we check if the condition (|9Tj) could be satisfied. Note Hq ~ Aj. Eq. ([95l) also indicates that in 
the inflationary era, where A = R ~ Aj , the magnitude of the mass is given by 



A2n+1 

2 



tin 
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which is large enough and the condition (|9"Tj) is satisfied. Here Eq. (|4"T)) is used. Thus, sufficient number 
of er-particles could be created. 

Let us consider the rotational curve of galaxy. As we will see the shift of the rotational curve does not 
occur due to correction to the Newton law between visible matter (baryon or intersteller gas) but due to 
invisible (dark) matter, and the Newton law itself is not modified. 

Let the temperature of the dark matter be T = \/k(3 where k is the Boltzmann constant. First, we 
assume the mass to ct of the scalar particle er is constant. As the total mass of dark matter is much larger 
than that of baryonic matter and radiation, we neglect the contributions from the baryonic matter and 
radiation just for simplicity. We now work in Newtonian approximation and the system is spherically 
symmetric. Let the gravitational potential, which can be formed by the sum of the dark matter particles, 
be V(r). Then the gravitational force is given by !F(r) = —mdV(r)/dr. If we denote the number density 
of the dark matter particles by n(r), in the Newtonian approximation, by putting re 2 = 8ttG, one gets 

2 pT 

T(r) = -^S- / Ans 2 n{ S )ds (98) 
r Jo 

and therefore V(r) is given by 



/T pS 
— J u 2 n(u)du 



(99) 



If one assumes the number density n(r) of dark matter particles could obey the Boltzmann distribution, 
we find 

n(r) = N o e- 0m " v(r) . (100) 

Here No is a constant, which can be determined by the normalization. Using (f99|) and (|100|) and deleting 
n(r), the differential equation follows: 

(r 2 V'(r))' = ATTGm a Nor 2 e~ Pm ' v{r) . (101) 
An exact solution of the above equation is given by 

V(r) = In (-) , r 2 = - 1 . (102) 

As the general solution for the non- linear differential equation (|101[) is not known, we assume V(r) could 
be given by (|102p . Then the rotational speed v of the stars in the galaxy could be determined by the 
balance of the gravitational force and the centrifugal force: 

V 2 7 77, 

m*— = -F{r) = m*V'{r) = . (103) 

r pm a r 

Here to* is the mass of a star. Hence, 

v 2 = — R , (104) 

that is, v becomes a constant, which could be consistent with the observation. 

For the dark matter particles from /(i?)-gravity, the mass m a depends on the scalar curvature or 
the value of the background a as in (|94p . The scalar curvature is determined by the energy density p 
(if pressure could be neglected as in usual baryonic matter and cold dark matter) and if we neglect the 
contribution from the baryonic matter, the energy density p is given by 

p{r) = m a n(r) . (105) 

Therefore it follows 

m a = m a (p(r)) = m a (m a n(r)} , (106) 
which could be solved with respect to m a : 

m a = rn a {n(r)) . (107) 
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Furthermore by combining (jlOOp and (|107p . one may solve m a with respect to V(r) and Nq as 

m a = m CT (JVb, V(r)) . (108) 

Then ([98]) could be modified as 

Hr) = ~ Gma (iV °' F(r)) f 47r« 2 m cr (JV„, F(r)) »(«)*, (109) 



which gives, instead of (|10ip . 

(rV(r))' = 47rGm ff (JV , K(r)) No ^ e -^(N ,v(r))v(r) (UQ) 



Eq. (|110p is rather complicated but at least numerically solvable. 

For the model (|46p . if the curvature is large enough even around the galaxy, the mass m a is given 
by (|95p . The scalar curvature A = R is proportional to the energy density (since the pressure could be 
neglected) , A oc p, and the energy density p is given by (|105[) . Then 



n|fi ^| 2(„ + 1 )(2„ + l)/„ |^ (iMr)ra . 



Using (llOOp , one also gets 



y(r) ,„^i), ^ s(Al) -ifi( wijft }"' (112) 

(2n + 3)/3m CT (r) m [ f{ J 

Here mo has mass dimension. By substituting (|112p into (|110p . it follows 



2n + 3/ /3 1 \ mo / m CT (r) 2 \ mo / m CT (r) 3 

+2r (l_ln!!kW'\£W| = lf 2(„ + l)(2n + !)/„ | A ^ 
V m / m CT (r) 2 J 2 { f( J 

It is very difficult to find the exact solution of (|113p . although one may solve (| 1 1 3[) numerically. Then we 
now consider the region where m„ -C mo but m(m CT /mo) is slow varying function of r, compared with 
the power of r. In the region, we may treat ln(m CT /mo) as a large negative constant: 

In K/ro ) ~ -C . (114) 

Then the following solution is obtained: 



m a (r) = m — 



2(2n + 3) 
2n+5 



r, 



2 



4(2n+l)(2n + 9)C , 2 ,|=±| f 2(n + l)(2n + l)/ 



(2n + 5)/3 v u; I fl 
Note that ro can be real for any positive n. Eq. pi2p shows that 

2(2n + 3) 

2(2n+ 1) 1 / r \ 2 ™+ 5 r 



= To - ^ - ■ (H6) 

2n + 5 pmo V r o/ 

Note that the potential (|116p is obtained by assuming the Newton law by summing up the Newton po- 
tentials coming from the /(i?)-dark matter particles ('scalaron') distributed around the galaxy. Eq. (|115p 
indicates that the condition m a <C mo requires r 3> ro. Then by using the equation for the balance of 
the gravitational force and the centrifugal force, as in (|103[) . we find 



(117) 
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which is monotonically increasing function of r and the behavior is different from that in (|104p . If there is 
only usual baryonic matter without any dark matter, the velocity is the decreasing function of r, if there 
is also usual dark matter, as shown in (|104|) , the velocity is almost constant, if dark matter originates from 
/(-R)-gravity, as we consider here, there is a region where the velocity could be an increasing function 
of r. Of course, one should be more careful as these are qualitative considerations. The condition 
m a <C mo requires r > but in the region faraway from galaxy, the scalar curvature becomes small 
and the approximation ()95|1 could be broken. Anyway if there appears a region where velocity is the 
increasing function of r, this might be a signal of /(i?)-gravity origin for dark matter. For more precise 
quantitative arguments, it is necessary to include the contribution from usual baryonic matter as well as 
to do numerical calculation. In any case, it seems very promising that composite particles from viable 
modified gravity which unifies inflation with late-time acceleration may play the role of dark matter. 

7 Discussion 

In summary, we reviewed i 7 '(i?)-gravity and demonstrated that some versions of such theory are viable 
gravitational candidates for unification of early-time inflation and late-time cosmic acceleration. It is 
explicitly shown that the known critical arguments against such theory do not work for those models. In 
other words, the modified gravity under consideration may pass the local tests (Newton law is respected, 
the very heavy positive mass for additional scalar degree of freedom is generated). The reconstruction 
of modified F(R) gravity is considered. It is demonstrated that such theory may be reconstructed for 
any given cosmology. Moreover, the partial reconstruction (at early universe) may be done for modified 
gravity which complies with local tests and dark energy bounds. This leads to some freedom in the choice 
of modified gravity for the unification of given inflationary era compatible with astrophysical bounds and 
dark energy epoch. As a final very promising result it is shown that modified gravity under consideration 
may qualitatively well describe dark matter, using the composite scalar particle from F(R) theory and 
Chameleon scenario. 

Thus, modified gravity remains viable cosmological theory which is realistic alternative to standard 
Einstein gravity. Moreover, it suggests the universal gravitational unification of inflation, cosmic accel- 
eration and dark matter without the need to introduce any exotic matter. Moreover, it remains enough 
freedom in the formulation of such theory which is very positive fact, having in mind, coming soon precise 
observational data. 
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